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ERRATA 

Page  10,  eqn  (25)  for  read 

2 

Page  20,  eqn  (38)  second  integral,  read  expi( -t  /2) 

Page  21,  eqn  (43)  for  k  read  K 

eqn  (47)  E(z^)  =  -34^j  4^^ 

eqn  (48)  for  4*2  read  4^^ 

Page  22,  eqn  (53)  E(z')^  =  -24^ 

line  I  The  skewness  is  not  zero,  and 
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Introduction 

The  present  state  of  theoretical  knowledge  of  wind  generated 
gravity  waves,  when  studied  as  a  random  process  in  nature,  is  based  on  a 
short  crested  Gaussian  sea  surface  obtained  by  linearizing  the  Eulerian 
equations  of  motion  (Pierson  [1952],  [1955];  Longuet -Higgins  [1957]),  on 
a  second  order  extension  of  the  Eulerian  equations  for  long  crested  waves 
in  which  the  underlying  linear  waves  are  Gaussian  (Tick  [1959]),  and  on 
various  results  involving  wave  spectra  that  do  not  specify  the  probability 
structure  of  the  waves,  some  of  which  involve  nonlinear  properties  of  the 
waves  (Phillips  [1958],  [1961]). 

The  short  crested  Gaussian  model  has  proved  fairly  useful 
in  explaining  the  angular  spreading  and  dispersion  of  swell,  the  marked 
variation  in  height  from  wave  to  wave  in  a  rtorm  sea,  wave  refraction, 
bottom  pressure  fluctuations  caused  by  waves  passing  overhead,  and  the 
motions  of  ships  in  waves  (St.  Denis  and  Pierson  [1953],  Lewis  [1955], 
Cartwright  and  Rydill  [1956]). 

These  applications  are  dependent  on  either  the  gross  features 
of  the  waves  or  on  natural  processes  that  appear  to  linearize  the  system 
even  further.  For  example,  the  wave  profiles  are  in  a  sense  averaged 
over  the  length  and  beam  of  a  vessel  at  sea,  and  the  nonlinear  effects  are 
reduced.  The  higher  harmonics  in  a  nonlinear  profile  and  the  high  fre¬ 
quency  linear  components  are  also  attenuated  with  depth  in  such  a  way 
that  the  application  of  Gaussian  noise  theory  to  the  zero  crossings  of 
bottom  pressure  fluctuations  leads  to  useful  results  (Ehrenfeld,  Good¬ 
man,  et  al  [1958] ). 

The  spectrum  of  the  time  process  at  a  fixed  point  has  been 
most  frequently  studied.  Since  the  actual  wave  is  recorded,  more  or  less, 
the  full  nonlinear  motion  is  recorded,  and  the  spectrum  of  this  nonlinear 
motion  is  estimated.  The  accuracy  of  computations  in  a  linear  theory 
based  on  spectra  estimated  from  a  nonlinear  record  is  open  to  question. 

In  particular,  higher  moments  of  the  spectra  are  quite  open  to  question. 

Theoretical  developments  in  the  study  of  the  short  crested 
Gaussian  sea  surface  require  numerous  higher  moments  of  the  linear  part 
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of  the  process.  If  computed  spectra  are  used  to  estimate  these  moments, 
many  of  the  higher  moments  are  obscured  by  high  frequency  white  noise  in 
the  estimate  or  made  questionable  by  the  high  frequency  nonlinear  part  of 
the  spectrum  so  that  the  results  predicted  from  them  are  doubtful.  If  the 
various  theoretical  forms  for  the  spectra  that  have  been  proposed  are  used, 
the  fourth  or  fifth  moment  of  the  frequency  spectrum  becomes  questionable 
and  the  meaning  of  the  second  or  third  moment  in  the  wave  number  spectrum 
is  obscure. 

These  problems  are  further  complicated  by  the  fact  that  sea 
surface  slopes  depend  on  the  very  high  frequency  capillary  waves.  With 
capillary  waves,  breakers,  whitecaps  and  turbulence  in  the  upper  layers  of 
the  water,  any  results  on  curvature,  as  they  depend  on  the  higher  moments 
of  the  spectrum,  are  highly  doubtful  for  a  storm  sea. 

The  actual  sea  surface  has  gross  features  that  are  not  in  accord 
with  the  sliort  crested  Gaussian  model.  No  artist,  for  example,  depicts 
wind  waves  as  irregular  sinusoidal  waves.  Wave  profiles  along  a  line  as  a 
function  of  distance  show  sharp  crests  and  long  shallow  troughs  that  do 
not  occur  in  the  model.  This  point  will  be  discussed  further  in  a  later  part 
of  this  paper. 

Purpose  of  paper 

The  purpose  of  this  paper  is  to  derive  three  new  random  sea¬ 
way  models  that  appear  to  have  some  properties  of  actual  seas  not  illus¬ 
trated  by  previous  models,  to  derive  a  few  of  these  properties  and  discuss 
other  possible  properties,  and  to  show  how  these  models  appear  to  be 
more  realistic  by  comparing  them  with  selected  observations.  These  new 
models  explain  some  of  the  difficulties  that  arise  with  the  higher  spectral 
moments.  They  may  also  be  capable  of  giving  some  information  about 
whitecaps,  breaking  waves,  and  sharp  crested  waves.  Since  all  mathe¬ 
matical  models  of  nature  fail  in  one  way  or  the  other  when  the  conditions 
and  assumptions  in  the  derivation  are  not  met,  these  models  will  provide 
a  second  choice  when  observations  of  waves  are  to  be  compared  with  theo¬ 
retical  results.  Some  questions  are  raised  that  can  only  be  answered  by 
additional  theoretical  work  and  by  precise  measurements  of  waves.  The 
theoretical  properties  of  these  models  can  be  compared  with  the  theoretical 
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properties  of  those  models  described  above  in  order  to  design  more 
definitive  experiments  concerning  the  nature  of  wind  generated  waves. 

The  Lagrangian  eqtiations 

Let  o,  p,  and  6  be  the  x,  y,  and  z  coordinates  of  a  particle 
of  fluid.  In  the  Lagrangian  system  of  aquations,  a  solution  to  the  equations 
consists  of  finding  the  positions  x,  y,  and  z  of  all  of  the  particles  in  the 
fluid  as  a  function  of  time  and  the  initial  positions  of  the  particles,  o,  p, 
and  6.  The  Lagrangian  equations,  according  to  Lamb  [1932]  are  given 
by  equations  (1)  where  subscripts  denote  partial  differentiation. 

^tt^a  +  ytt^a  +  <^t  ^  +  Pa^P  =  ° 

(1)  ^t^P  +  ^tt^p  +  <^t  "  g^^P  ^  Pp/P  = 

^tt^6  +  ytty6  +  <^t^g>^6  --Ps/p^o 

The  equation  of  continuity  is  most  conveniently  expressed  by 

(2)  d 

at 

Such  solutions  need  not  be  irrotational,  but,  if  a  function, 

F(a,  p,  6,  t)  ,  can  be  found  such  that 

(3)  dF  =  y^y^+z^zjda+  (x^Xp+  y^  yp  +  z^  Zp)dp  +  (x^  x^  +  y^  y^  +  (z^  Z5)d6 

is  a  perfect  differential,  there  is  no  vorticity. 

A  zero  order  solution  to  these  equations  is  given  by 

X  =  a 

y  =  P 

z  =  6 

p  =p^  -  g  p  6 

iri  which  all  fluid  particles  are  at  rest  in  hydrostatic  equilibrium  under  the 
force  of  gravity. 

We  expand,  following  the  concepts  described  by  Stoker  [1957]  . 
about  the  zero  order  solution  in  terms  of  a  small  parameter,  e  ,  as  in 
(5)  in  which  F^  is  a  constant.  Here  we  think  of  t  as  equal  to  ak  , 


9(x,  y,  z) 
8(a.p.6), 
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but  the  first  order  terms  must  have  the  dimensions  of  a  length,  and  so 
a  =  t/k  is  used  in  the  various  solutions  that  are  obtained.  The  parameter 
t  never  appears  explicitly.  See,  for  example,  Longuet -Higgins  [1953], 
and  Pierson  and  Fife  [1961]  for  different  ways  that  e  can  be  interpreted. 

X  X  2 

X  =  a  +  eXj  +  e 

,5,  ^  '  P  *  ‘V,  *  .^2 

z  =  6  +  £Zj  +  e^Z2 

P=P^-gp6+EPj  +£^2 

F  =  F  +  £F,  +  £^F, 
o  1  2 

The  equations  in  e  become 

^Itt^  Piq/P  =  ° 

yitt^8^ip^Pip/p  =  o 
(6)  ^Itt^8^16+Pl6/P  =  0 


^lat  +  yipt+^16t  =  ° 

dFi  =  Xj^da  +  yj^dp  +  Zj^d6 
2 

The  equations  in  e  become 


(7) 


^2tt  ^  8^2a  ^  PzJP  =  -  ^in^la  '  yitt^la  '  *ltt  "la 
y2tt  '''  ®"2p  P2p'^P  "  ■  ^Itt’^ip  ■  ^Itt^ip  ■  ‘"ltt"ip 
"^2tt  ®"26  P26^P  ""  ■  ’‘itt’^ie  ■  7ltt7’l6  '  "ltt"l6 

f^2a^y2p^"26  +  yip^la''’^la"l6^yip"l6-"loyi6 
-"ipyia-^ipyiaJt  =  ° 

^<"2t^"ltyi6''yityi6'-"lt"l6>‘*® 
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Equations  (6)  are  linear.  A  solution  of  them  in  turn  determines  the 
right  hand  side  of  (7).  Since  the  left  hand  side  of  (7)  is  linear,  in  principle 
at  least,  a  solution  can  be  found.  Also  in  principle  (though  not  in  practice), 
it  is  possible  to  proceed  to  as  high  an  order  as  desired  by  this  procedure. 


The  Gerstner  wave 

A  solution  to  equations  (6),  applicable  to  waves  in  deep  water,  when 
added  to  the  zero  order  solution,  yields  equation  (8). 

X  s  a  +  exj  =  a  -  sin(kja  +  -  wt) 

^■^7  If  fi 

=  P  — ic  8in(k^a  +  -  t*Jt) 

(8)  z  =  6-(-£2j=6  +  a  e^^  co8(kjO  -  k2P  -  wt) 


P  =  Pq  +  ePl  =Po  -  gp6 

Fi  =  sin(kjO  +  k^P  -  wt) 

In  order  to  impose  the  condition  that  p(o,  p,  6)  =  p^  at  6  =  0,  it 
is  required  that 


(9) 


k 


The  continuity  equation  imposes  the  condition  that 


(10)  k^=kj^+k2^ 

The  first  order  solution  has  no  vorticity  to  first  order  as  F ^  has 
been  found. 

An  alternate  form  for  the  solution  given  by  (8)  is  (11). 

2*/  2 

X  =  a  -  a  cos  9e“  '  ®  sin(  —  (a  cos  0  +  p  sin  0)  -  ut) 

S 

2c/  2 

y  =  P  -  a  sin  0e“  '  ®  0  +  P  sin  0)  -  «t) 

2./  2 

z  =  6  +  a  e“  cos(^  (o  cos  0  +  p  sin  0)  -wt) 


(II) 
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These  solutions  (eqs.  8  and  11)  are  simply  disguised  versions  of 
the  Gerstner  wave  in  deep  water  as  described  in  Lamb.  The  angle  0  is 
the  direction  toward  which  the  wave  is  traveling.  When  the  solution  is 
periodic,  as  above,  the  original  nonlinear  equations  are  satisfied  exactly, 
but  it  can  be  shown  that  there  is  some  second  order  vorticity  and  that 
there  is  a  second  order  correction  to  the  pressure  given  by 


/loi  e  .  a^k  2k6  a^k 

(12)  P  =  Pq  ■  g  P  0  +  g  P  *  ’“j^gP 

Another  possibility  also  suggests  itself.  With  =  a^ukt  e^^^  cos  0, 

y^  =  a^ukt  e^^^  sin  0  ,  and  =0,  equations  (1),  (2),  and  (3)  are  still 
satisfied  to  second  order  but  third  order  complications  arise.  These  two 
terms  are  the  equivalent  of  the  s  econd  order  current  found  in  irrotational 
waves  in  the  Eulerian  system  of  equations  as  described  in  Lamb. 

Extension  of  these  results  to  finite  depth  would  probably  follow  the 
results  of  Biesel  [1952]  who  achieved  a  remarkable  representation  for 
periodic  breaking  waves  by  means  of  the  Lagrangian  equations  (see  also 
Pierson  [1955]).* 


A  randomized  short  crested  model 

The  expressions  for  x^.y^.and  z^  in  either  (8)  or  (11)  are  solutions 
of  equations  (6).  These  equations  are  linear.  Therefore,  a  sum,  over  dif¬ 
ferent  parameter  values,  of  solutions  of  the  form  of  x^y^,  and  z^  will  also 
be  a  solution.  If  it  is  assumed  that  each  pairticle  has  a  displacement  from 
its  rest  position,  described  by  a  stationary  Gaussian 

vector  process  with  specified  coherency  relationships,  a  solution  to  equations 
(1)  that  is  correct  to  first  order  is  given  by  equations  (13). 

«>  IT  2 , .  2  ___________ 

x  =  a-/  ®  cos0sin(— (aco80+  p  sin0(-wt  +  e(w,0))  n/2S(w,0)  dOdw 

o  __  g 


TT  2.  2 

(13)  y  =  p-7  /e"  ^^®8in0  8in(—(aco8  0+ Paine)  -«t  +  «(w,0)) '«i2S(u.0) d0d« 
o  -v  * 

If  2  2 

z  =  6+7  /  e***  cos  (  —  (aco80+ p  8in0) -«t  +  c(w,0))  '^2S(u,0)  d0du 

O  -IT  ® 


*rhe  remarks  concerning  the  Gerstner  wave  are  in  error  as  these  results 
six  years  later  show. 


-7- 


The  notation  here  ia  that  need  by  Pierson  (1952, 1955)  except 
S{u,  6)  is  the  resolution  of  the  variance  of  the  particle  motions  into  fve* 
quency  and  direction.  For  tome  fixed  particle,  say,  one  at  the  suriface 
given  by  a  s  0,  p  =0,  6=0,  the  motions  x  =  x(t),  y  =  y(t),  and  z  =  z(t) 
form  a  stationary  vector  Gaussian  process  such  that  the  spectriun  o£ 
x(t)  is 

/s(a).  e)(cose)^de , 
the  spectrum  of  y(t)  is 

/$(u,  e)(sine)^de , 
the  spectrum  of  z(c)  is 

/SC«.  6)  d0  . 

the  cospectrum  of  x(t)  y(t)  is 

/S(u,  6)  cos  6  sinO  d6  , 
the  quadrature  spectrum  of  x(t)  z(t)  is 
/S(u,  6)  cos  6  dO  , 

the  quadrature  spectrum  of  y(t)  z(t)  is 
/S(u,  6)  sine  de  , 


and  all  other  cross  spectra  are  zero. 

Other  notations  are  often  useful.  For  example,  a  notation  sinxilar 
to  that  of  Longuet -Higgins  [1957]  would  yield  (14)  as  the  complete  equi¬ 
valent  of  (13). 


(14) 


_  m  k6 

X  =  o  -  £  a  ,  -r— e 
mri 


y  =  p-2a 


sin(k^ 

8in(k^ 


a+kp-w  t+€ 
n  mn 


a  +  kp-a>  t  +  € 
rT  mn 


mn 


mn 


) 

) 


(15) 


z  =  6  +  r  a  e'^®  co8(k  a  +  k  w  -  w  t  +  «  ) 

mn  m  n  mn  mn 

In  (14),  added  conditions  are  that 


=  k 


2 

mn 


and  that 
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(i6i 


2 


“rnn  _  XT? 

g  "V  m 


The  free  surface 

Those  particles  on  the  free  surface  are  described  by  the  condition 
S  =  0  and  (13)  becomes 


J  /sin(— (acos6  +  p  sln0)-wt  +  «)  cos6  JZS{u,  8)  d6du 
o  -ff  * 

(17)  y  =  y(o,p,t)  =  p  -  /  /sin(^i-(acos0  +  p  sin©)- wt  +  «)  sin6  V2S(w,  0)  d6du 

O  -IT  * 

«0  IT  2 

z  =  z(a,p,t)  =  /  /  cos(—  (a  cos0  +  p  sin0)  -  wt  +  *)  •/2S(w,  0)  d0  dw 
o-ir  ® 


In  principle,  and  perhaps  admitting  triple  values  for  the  inverses 
over  limited  regions  in  the  x,  y,  t  space,  x  =  x(a,  p,  t)  and  y  =  y(o,  p,  t) 
imply  inverses  of  the  form 


(18) 


a  s  a(x,  y,  t)  and 
p  =  P(x,  y.  t) 
so  that  the  free  surface  can  be  found  from 
(19)  z  =  z(o(x,  y,  t),  P(x,  y,  t)t)  =z(x,  y,t) 


The  surface  so  defined  is  certainly  not  the  equivalent  of  the  short 
crested  Gaussian  sea  surface.  The  short  crested  Gaussian  sea  surface 
is  equivalent  to  this  representation  when  the  amplitude  of  the  particle 
motions  becomes  so  small  that  a  =  x  and  P  =  y  are  satisfactory  approxi> 
mations  to  the  inverses  given  by  (18).  Otherwise,  further  study  suggests  that 
z  =  z(x,  y,  t)  has  many  features  of  actual  waves  in  that  the  crests  can  be 
quite  sharply  pointed  and  in  that  the  higher  nonlinear  harmonics  that  occur 
in  the  higher  order  derivations  in  the  Euler ian  system  are  already  present 
in  this  model. 

It  is  believed  that  the  problems  that  arise  in  the  adequate  probabi¬ 
listic  description  of  the  surface  defined  by  (19)  in  terms  of  (17)  and  (18) 
will  be  very  difficult  to  solve  and  that  considerable  effort  will  be  required. 
All  of  the  problems  that  arise  in  the  study  of  Gaussian  noise  and  in  the 
study  of  the  short  crested  Gaussian  sea  surface  have  their  analogue  in  this 
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model. 

Second  order  effecta  in  the  abort  cre>ted  model 

For  reasons  to  be  discusaed  later,  it  may  be  necessary  to  proceed 
to  second  order  in  this  model  and  determine  the  effects  of  equations  (7)  on 
the  wavy  surface  especially  if  breaking  waves  are  to  be  studied. 

The  long  crested  linear  model 

If  S(u,  6)  becomes  concentrated  at  the  angle  6  =  0  and  degenerates 
to  a  function  of  u  only,  equations  (13)  become 

2  2 

X  =  x(a,  6,  t)  =  a  -  /  e*^  8in{—  a  -  cat  +  f)  ^/2S(ca)  dZT 

o  8 

(20)  y  =  y(P)  =  p 

00  2-,  2 

z  =  z(o,  6,  t)  =  6  +  /  e*^  '  ®  co8(—  a  -  cot  +  €)  \lZS{u)  dca 
o  * 

An  alternate  notation  i8  given  by 

k  6 

X  =  a  -  £  a  e  8in(k  a  -  ca  t  +  c  ) 
n  n  n  n 


z  =  6-  2:a  e”  co8(k  o  -  w  t  +  c  ) 
n  n  n  n 

For  6  equal  to  zero,  the  inverse  of  the  first  equation  is 

(22)  a  =  a(x,  t) 

This  implies  a  free  surface  given  by 

(23)  z  =  z(o(x,  t),t)  =  z(x,  t) 

Second  order  effects  in  the  long  crested  model 

The  long  crested  model  is  more  amenable  to  an  investigation  of 
second  order  effects  because  it  is  simpler.  We  consider  the  problem  of 
two  waves,  and  then  generalize  to  the  randomized  process.  For  two  waves 
the  linear  solution  is  given  by 


iFgaaaiMIMtjw 
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kjO  k,6 

Xj  =  -aje  *  sin(k^a  -  u^a  +  <j)  -  a2e  ^  •in(k2a  -  <*>2*  +  <3^ 

k.6  k26 

(24)  *1  “  ®’i  ®  ^  co8(k^a  -  Wjt  +  «j)  +  a2e  cos(k2a  -  ci)2t  +  C2) 

Pj  =  0 

a.u,  k,8  ^9^ 

Fj  =  "  e  8in(k^a  -  Wjt  +  6  8in(k2a'-  ^2) 


It  i8  irrotational  to  fir8t  order.  For  8implicity,  we  a88ume  that  U2  >(o^. 

The  eecond  order  eqviatione  to  be  eolved  are  obtained  by  subatitut- 
ing  (24)  into  (7)  where  y  is  set  equal  to  ^  only  and  all  partials  with 
respect  to  p  vanish.  These  equations  are  given  by 


*2tt  +  8*2a  ^  P2a  ^ 


2k,  6  ^  o  2k2S 


*2t  ^  «*26  P26/P  =  ^  +  ^2%^  * 

(kj+k2)6 


+  2aja2kjk2e 


cos((k2-kj)a-(w2-Wj)t  +  «2-«j) ) 


(k,+  k2)6 

*26  coe((k2-kj)o-(w2-Wj)t+ €2-«i)Jt  *  0 


-  2k.(».  2k,6  ,  ^^9^ 

dF2  =  [X2j  -  (aj  kjW^e  ^e  +  e 

(kj+k2)6 

+  aja2(Wjk2  +«2^)  ®  co8((k2  -kj)o-  (w2-Wj)t  +  C2  -€j) )  Jdo 

a.a2  (k^+k^lS 

+  [(*2t  ■  ^“g”  ^“2  ■  "l^^"2"l^  ®  8in((k2  -  kj)o -  -  «j)t  +  C2  -  «j)Id6 

In  (25),  we  have  four  equations  in  p,  z  and  x.  The  first  three  de¬ 
termine  a  solution,  subject  to  appropriate  boundary  conditions,  for  p,  z, 
and  X.  The  last  is  then  used  to  see  if  the  solution  is  irrotational  and  to 
make  the  solution  irrotational  if  possible.  Note  that  x  =  X2(6)t  can  be 
added  to  (23),  if  desired,  as  part  of  the  total  solution  to  the  equations.  The 
boundary  conditions  are  that  P2  =  0  at  6  =  0  and  that  X2  and  <2  approach 
zero  as  6  approaches  -  <*  . 

If  the  right  hand  side  of  the  second  equation  is  represented  by 
g[F(6)  +  G(a,  6,t)]  and  the  third  equation  is  written  as  X2^'*’ 82^ -G(a,  6,  t)  s  0, 
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it  ia  possible  by  cross  differentiating  the  first  two  equations,  subtracting 
and  combining  the  result  with  the  third  equation  to  obtain  the  result  that 

(26) 


*2ttaa  *2tt66  “  '^tt6 


=  +  8G, 


66 


This  provides  the  inhomogeneous  solution  for  Z2(a,  6,t).  The 
third  equation  then  determines  X2(a,  6,t)  and  the  first  two  yield  the  same 
solution  for  P2- 

Howe  ^er  P2  does  not  satisfy  the  free  surface  condition,  and  a 
solution  to  the  homogeneous  equations  must  then  be  found  that  when  added 
to  this  solution  yields  a  proper  form  for  P2> 

This  provides  solutions  to  the  first  three  equations.  When  ^2 
Z2  are  substituted  into  the  last  equation,  it  is  seen  that  the  vorticity  can 


2  2kj6  2  2^2* 

be  made  zero  by  adding  a^  kjtOjC  t  +  a2'^k2W2*  *  *2‘ 


The  second  order  terms  therefore  become 


3  3 

a, a,  w,  +w,  \  (k-+k,)6 

'  je  sin((k2-kj)tt *2'*1^ 


/“l  ^^2  \ 

V  W2-“/ 


^1^2  {1^2  "^1)^ 

+  -j-=-  (W2  +  e  sin((k2  -  kj)a  -  *2  '  *  1  ^ 

2  ,  2  ,  ^^2^ 

+  aj  Wjhj  e  f  *2  “2^2  ®  * 


^1^2  2 


(27)  z,  =  — — (wj  +  «jW2  +  «2 


2  (kj+k2)6 


g 


coe((k2-kj)o  -  (w2  -  tOj)t  +  €2  -  c^) 


{k2-kj)6 

^  «2(«2'''‘^)«  cod(k2-kj)a-(w2-«j)t  +  €2  -  Sj^) 


*i»2 


a,^k,  2k,  6  a-^k,  2k- 6 

P2=gp4^(e  ^  -l)  +  gP-V^(e  -1) 

(k.+k2)6 

-  2paja2W2“i  «  cod(k2-kj)o-{«2‘"i)*  ^  *2"*l) 

{k2-k.)6 

+  2paja2W2“i  «  coi((k2-kj)o- {«2<*>iH  *2‘‘i) 

F2  =  a^a2W2®^^^'*'^^*^“^^^2"  V“”^"2"“l^*  ‘2’‘l^ 

(k2+kj) 


+  a2aj(.^-a.j). 


siiiCkg  -  kj)o-{«2-«iH  +  *2"*l«^ 


The  full  lolution  ie  obtained  by  combining  equations  (5),  (24), 
and  (27).  This  solution  satisfies  the  equations  to  second  order.  It  is 
irrotational  to  second  order. 

If  more  terms  are  added  to  the  linear  solution  subject  to  the  con¬ 
dition  that  u,  <  u~  <  u.  <  •  ‘  •  <  tt)  ,  the  terms  in  the  linear  solution  inter- 
i  £  i  n 

act  in  a  predictable  way  to  generate  appropriate  second  order  terms. 

The  randomized  second  order  solution  for  x(a,  6,t)  and  z(a.  6,t)  are 
given  by  equations  (28). 


x(a,6,t)  =  o  -  £a^e  sin(k^a-u^t  + 

/  (k,+k.)6 


j>ii  e  \  “j 

a.  a  (k.-k.)6 

+  £  Z - «(w.  +  w.)w,e  ^  8in((k.-k.)o-(w,-w,)t  +  €.-€.) 


(28) 


2  2k.  6 

+  Z  a.^u.k.  e  ^  t 
111 

k,6 


z(a,6,t)  =  6  +  Z  a.  e  ^  cos(k^a-b>^t  +.  €^) 

a.  a,  2  2 

+  Z  Z-^  (w.  +«.Wj  +  «j  )  e  ^  cos((kj-kj)o-(Wj-Wj)t  +  fj-e  j 


J>i 

Z 

J>i 


a.  a.  (k.-k.)6 

-  Z  Z Wj(«j  +  Wj)  e  J  cos((kj-kj)o-(<.>j-<Oj)t+ €j-€.) 


The  solution  given  by  equation  (28)  has  features  that  are  quite 
different  from  solutions  obtained  by  the  assumption  of  irrotationality  in 
the  Zulerian  system.  The  trigonometric  second  order  terms  involve  only 
the  difference  between  two  frequencies.  One  of  the  second  order  terms 
dies  out  rapidly  with  depth.  The  other  dies  out  slowly  with  depth.  The 
term  that  is  linear  in  t  is  the  average  drift  of  the  particle,  and  as  a  group 
of  higher  waves  passes,  the  effect  of  the  second  order  terms  is  to  increase 
the  drift.  Correspondingly,  as  low  waves  pass,  the  drift  is  decreased. 

The  positions  of  those  particles  that  are  on  the  free  surface  are 
obtained  by  setting  5=0,  and  the  result  is  equation  (29). 


(29) 


x(o,t)  =  a  -  r  aj  •in(k^a  -  Wj 


-r  I 
j>i  i 


“j-“i 


t  +  «i) 

2 

■  iin((kj-  kj)a  -  (u^  -  «»j)t  +  €j  -  f^)  +  X  a^  w^kjt 


*(o,t)  =  T  a.  coa(k^a>(i)^t  +  f.) 

+  _X  r  !^a.2coa«y  kj)a-(a)j.a>j)t  + 


The  parametric  equations,  x  =  x(a, t)  and  z  =  s(a,  t),  imply 
an  inverse  for  x  =  x(a.t)  such  that  a  =  a(x.  t),  and  so 

(30)  *  =  a(o{x,t).t) 

is  the  equation  for  the  free  surface  correct  to  second  order  as  obtained  by 
this  solution  of  the  Lagrangian  equations  of  motion. 

Comparison  of  the  Lagrangian  and  Eulerian  models 

Results  obtained  in  the  Lagrangian  system  of  equations  are 
difficult  to  compare  with  results  obtained  in  the  Eulerian  system  of 
equations.  In  the  Eulerian  system  the  fluid  velocity  at  a  fixed  point  below 
the  surface  will  be  Gaussian  in  the  linear  Gaussian  model.  The  velocity 
at  a  fixed  point  in  the  linear  Gaussian  model  in  the  Lagrangian  system  has 
not  been  found  but  it  involves  finding 

dx 

U  =  gj- =  u(o(x,  Z.t),  6(x,  Z,t),  t), 

(31)  a  =  a(x,  z,t)  ,  and 
6  =  5(x,  z,t) 

from  equations  (21),  even  in  the  long  crested  case,  and  this  velocity  is 
not  Gaussian  unless  the  waves  are  so  low  that  a  =  x  and  6  =  z  are  suf¬ 
ficiently  accurate  inverses  of  x  =  x(a,  6,  t)  and  z  =  z(a,  6,t). 

However,  a  higher  order  model  in  the  Eulerian  system  of  equations 
may  well  yield  velocities  comparable  to  those  given  by  a  model  to  a  lower 
order  in  the  Lagrangian  equations. 

The  first  and  second  order  models  in  the  Lagrangian  equations 
appear  to  have  advantages  over  the  comparable  first  and  second  order 
models  in  the  Eulerian  equations,  as  will  be  shown  later.  Since  the 
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resulting  free  surface  is  non-Gaussian  it  may  be  Just  as  difficult  to  obtain 
results  on  the  probability  structure  of  the  first  and  second  order  Lagrang- 
ian  models  as  it  would  be  to  obtain  results  on  a  much  higher  order  Eulerian 
model.  If  the  problems  prove  intractable,  numerical  computation  is  a 
possibility  that  should  yield  some  results. 

The  real  problem  is  thus  to  compare  the  six  available  random 
models  with  nature  and  to  devise  relevant  experiments  to  see  which  model 
comes  the  closest  to  agreeing  with  that  which  is  observed.  Any  linear  non- 
Gaussian  model  that  can  be  devised  along  with  its  higher  order  extensions 
should  also  be  allowed  to  enter  the  competition.  We  suspect  that  such 
models  will  not  do  as  well. 

Re  c  apitulation 

So  far  in  this  paper,  three  new  random  sea  surfaces  have  been  obtain¬ 
ed.  They  are  the  short  crested  sea  surface  given  by  (17),  (18),  and  (19) 
based  on  a  linear  superposition  of  the  particle  motions;  the  long  crested 
sea  surface  given  by  (21),  (22).  and  (23)  based  on  a  linear  superposition  of 
the  particle  motions;  and  the  long  crested  sea  surface  given  by  (29)  and 
(30)  based  on  the  second  order  correction  to  the  long  crested  model.  Each 
model  is  irrotational  to  the  order  to  which  it  is  carried  out.  The  notation 
chosen  to  represent  the  second  order  long  crested  model  is  not  as  useful 
as  the  notation  used  by  Tick  [1959]  in  his  study  of  the  second  order  long 
crested  model  in  the  Eulerian  system.  These  results,  however,  show  the 
existence  of  such  a  model. 

Properties  of  the  models 

The  multivariate  probability  structures  of  these  models  have 
not  yet  been  found.  Only  a  few  of  the  properties  of  the  long  crested 
linear  model  have  been  found.  A  possible  realization  of  z  =  z(x)  has 
been  constructed,  and  the  probability  density  of  z(x)  for  a  fixed  x  has 
been  found. 
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Comtruction  of  g  =  z(x)  for  the  linear  model 

With  6  and  t  equal  to  zero  in  (21) ,  the  result  is 

x  =  a+x.  =  a-  Xa  sin(k  o  +  t  ) 

(32)  ^  n  n  n 

z  -  z,  =  Za  cos(k  a  +  e  ) 

1  n  n  n 

The  quantities  x^  and  z^  form  a  vector  Gaussian  process  *  The  spectrum 
of  Xj  is  S(k),  the  spectrum  of  z^  is  S(k).  the  co -spectrum  is  zero,  and 
the  quadrature  spectrum  is  -S(k).  The  coherency  is  one. 

Tick  [19(>1]  has  developed  procedures  for  the  construction  of  vector 
Gaussian  processes  with  required  spectra,  cross  spectra,  and  coherencies. 

A  vector  process  with  a  reasonable  spectrum  and  with  the  proper  cross 
spectra  was  available,  and  it  was  a  simple  matter  to  prepare  graphs  of 
x(a)  and  z(a)  as  a  function  of  a. 

The  results  are  shown  in  figure  1,  where  x^{a),  z(a),  and  al-x^(a) 
are  all  graphed  as  a  function  of  a.  For  any  value  of  a  in  this  figure,  a 
pair  of  values  for  x  and  z  result  that  can  be  plotted  as  a  point  in  the  x,  z 
plane  as  shown  immediately  below.  The  points  in  the  x,  z  plane  then  trace 
out  the  indicated  curve  that  represents  z  =  z(x).  Only  two  crests  are  shown 
in  figure  1. 

Figure  2  shows  the  result  of  preparing  much  longer  graphs  of 
x(a)  and  z(a).  Eight  crests  are  shown  in  the  figure.  In  many  ways,  this 
artificial  record  appears  to  be  much  more  realistic  than  one  that  might  be 
obtained  with  a  linear  Gaussian  model. 

These  results  show  that,  although  the  spectrum  of  the  orbital  motion 
could  be  band  limited  (i.  e.  ,  identically  zero  outside  of  a  certain  wave  num¬ 
ber  range),  the  spectrum  of  z  =  z(x)  could  be  very  rich  in  higher  harmonics 
of  the  orbital  motion  spectrum.  With  one  or  more  simple  discontinuities  in 
slope,  as  would  be  achieved  at  the  sharp  crests,  the  asymptotic  form  pre¬ 
dicted  by  Phillips  [1958]  would  be  found  as  applied  to  the  long  crested  case. 
However,  it  would  occur  at  wave  numbers  much  higher  than  are  presently 
capable  of  being  resolved  in  spectral  analyses. 

One  of  the  difficulties  with  the  Eulerian  models  is  that  they  do  not 
indicate  within  themselves  unrealistic  wave  forms.  A  spectrum  that  would 
yield  an  impossible  sea  condition  does  not  appear  to  be  different  from  one 
that  would  yield  a  realistic  sea  condition.  The  linear  Lagrangian  model  is 

*In  fact,  Zj  is  the  Hilbert  transform  of  x^. 
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1.  GRAPHS  OF  Z(a),X,(a),  AND  Xia)  TO  SHOW  HOW  Z  •  ZOO  IS  OBTAINED 
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already  capable  of  differentiating  between  completely  unrealistic  apectra 
and  more  realistic  ones.  If  the  space  scale  of  and  is  changed  for 
the  same  vector  realisation  used  to  construct  figures  1  and  2,  the  result 
is  figure  3.  Here  o  is  triple  valued  for  a  certain  range  of  x,  and  when 
s(x)  is  plotted  the  surface  forms  a  large  loop.  This  may  be  indicative  of  a 
very  unstable  state  in  the  waves  that  could  not  be  attained  in  nature  be* 
cause  breaking  waves  would  form  each  time  such  a  loop  would  be  initiated 
in  the  space -time  representation  and  the  kinetic  energy  of  the  wave  motion 
would  be  dissipated  by  turbulence  in  the  whitecaps  to  an  extent  that  the 
particle  motion  spectra  would  be  modified  to  a  more  realistic  condition. 

Conditions  for  a  loop  in  z  =  z(x) 

A  loop  occurs  in  z  =  z(x)  if  z  is  triple  valued  as  a  function  of  x. 
This  in  turn  occurs  if  x(a)  is  not  monotonically  increasing.  Thus  if 
m  <  0  where 


(33)  m  =  — *  1  -ra^k^cos(k^o.w„t+€^) 

a  loop  will  form. 

The  probability  that 

(34)  £a  k  co8(k  o-wt  +  t  )>1 

n  n  n  n  n 

must  therefore  be  investigated. 

The  three  moments  given  by 


oc 


(35) 

vP  =  /  S(k)  dk 
o  o 

(36) 

=  /S(k)  k  dk 

and 

(37) 

=  7  S(k)k^dk 

^  O 

need  to  be  considered  in  this  connection,  and  from  (34)  it  is  evident  that 
must  exist  and  be  bounded. 

The  probability  that  m  will  be  less  than  zero  is  given  by 
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(38^  p(m<0)  a 


r 


dti 


J. 


■/2«1», 


dt 


Equation  (38)  may  be  interpreted  to  be  the  probability  that  aome  particle, 
a^t  will  be  found  at  the  top  of  a  loop  between  the  points  where  z  *  z(x)  iz 
vertical  for  some  particular  time  of  observation.  Such  a  particle  is 
surely  involved  in  a  breaking  wave.  Other  particles  that  do  not  satisfy 
this  condition  also  may  be  involved  in  a  breaker. 

Since  the  spectrum  S(k)  has  not  as  yet  been  measured,  it  is 
necessary  to  investigate  possible  relationships  between  S(k)  and  spectra 
of  z(x,  0)  and  z(0,t)  before  estimating  the  probability  in  (38). 

Marginal  distribution  of  zaz(x) 

The  graphical  construction  used  to  illustrate  the  generation  of  the 
random  function  z  =  z(x)  can  be  used  to  obtain  the  marginal  distribution  of 
z  =  z(x).  The  joint  density  of  z(a,  0)  and  ^  =  m(a,  0)  is  bivariate  normal. 

The  variance  of  z(a)  is  41  >  the  variance  of  m(x)  is  4,.  and  the  covariance 
is  as  defined  in  (34),  (35)  and  (36). 

If  f(z,  m)  is  the  bivariate  normal  distribution  of  z  and  m,  we 

note  that 


(39)  F(z,  m)  =  f(z,  m)m 

nearly  has  the  required  properties  of  a  probability  density  function  since 
m  is  dimensionless  and  £(m)  =  1.  Also  since  x  and  a  have  the  same 
dimensions,  a  given  number  of  equally  spaced  points  in  a  small  interval 
of  a  given  length  over  a  produce  m  times  this  number  of  points  with  this 
same  spacing  in  an  interval  of  a  different  size  on  the  x  axis. 

The  function  F(z,  m)  is  not  a  probability  density  function  strictly 
because  it  is  negative  for  m  <  0,  but  this  corresponds  physically  to  im¬ 
possible  configurations  of  z  s  z(x),  and  points  on  z  =  z(x)  generated  when 
m  <  0  will  be  shown  to  have  extremely  low  probabilities.  Actually 
even  more  points  on  z  =  z(x)  are  probably  destroyed  by  the  breaking  pro¬ 
cess,  and  z  =  z(x)  is  modified  in  form  long  before  m  becomes  zero  for 
some  sea  conditions. 

However,  as  an  approximation,  we  can  truncate  F(s,  m)  at 
m  s  0  and  renormalize.  Thus 


(40) 
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and 

(51) 

(52) 

(53) 

(54) 


E(2')  >  0 
E(2'^)  = 

E(*'^)  =  0 

E(»'^)  =  34*^^  -  64<p4<j^  -  34' 


The  skewness  is  zero  and  the  coefficient  of  excess  is 
■64'iV(4'o-4'i2)^  • 

Estimates  of  *^2 

From  (52)  the  average  of  z'(x)  is  zero  and  the  variance  is, 

2 

say,  4'  *  =  4*  -  4',  •  This  is  the  variance  that  is  estimated  by  a  free  sur 
o  o  i  ' 

face  wave  record. 

Since 

(55)  S(k)  dk  =  S{^)  =  S«(w)  du 

9  9 

S  S"(f)df  . 


and  since  (for  example) 

(56)  S(k)  k^  dk  =  S"(f)  I6w^  df 

g 


let  us  assume  that 


(57) 


s"(f)  =1 


K/f”  f,  <  f  <  f 
I  u 

0  otherwise 


From  equations  (35),  (36),  (37),  (55),  (56),  and  (57),  one  can 
obtain  the  result  that 


(58) 

(59) 


g  o  l,  . 


n-3 

T 


and 
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(60) 

where 

(61) 


0=  (1. 


64t 

% 


w  *ri- 

■  Vt 

ro  Lj 

/n  - 

'r!  ■  'V'u'”"’] 

(rrrj 

Li  - 

♦  2 

For  one  fairly  severe  storm  at  sea  4^^  was  about  equal  to  116  ft 
and  fj  was  equal  to  0.046.  The  results  above  permit  an  investigation  of  the 
sensitivity  of  the  particle  motion  moments  to  the  values  of  n  and  f^. 

Tables  1  through  5  give  the  values  of  'I'j,  4'2>  1/*'/^^  » 

2 

*1',  /4'  4',  for  various  values  of  n  and  for  f  =  0.46,  0.92,  4.6  and  *  . 

1  o  c  u 

The  values  of  f^  correspond  to  periods  of  2.16,  1.08,  .216  and  zero  seconds. 


Table  1. 


Values  of  tl'  for  various  values  of  n  and  f 

O  ^  a*!  ^ 

for  4*  *  =  ll6ft^  and  f-  =  0.046  sec  . 
o  i 


/ 

'  c 

0.46 

0.92 

4.6 

00 

5.9 

116.26 

116.26 

116.26 

116.26 

5.7 

116.28 

116.28 

116.28 

116.28 

5.5 

116.30 

116.30 

116.30 

116.30 

5.3 

116.32 

116.32 

116.32 

116.32 

5.1 

116.35 

116.35 

116.35 

116.35 

4.9 

116.38 

116.39 

116.39 

116.39 

4.7 

116.42 

116.43 

116.44 

116.44 

4.5 

116.47 

116.49 

116.50 

116.50 

-24 


Table  2. 

Values  of  for  various  values  of  n  and  f^ 

for  'I'  *  *  116  ft^  and  f^  ■  0.046  sec”^. 
o  i 


f 

u 

n 

0.046 

0.92 

4.6 

00 

5.9 

-«s51— 

0.51 

0.51 

0.51 

5.7 

0.53 

0.53 

0.53 

0.53 

5.5 

0.55 

0.55 

0.55 

0.55 

5.3 

0.5? 

0.57 

0.57 

0.57 

5.1 

0.59 

0.59 

0.59 

0.59 

4.9 

0.62 

0.62 

0.62 

0.62 

4.7 

0.65 

0.66 

0.66 

0.66 

4.5 

0.69 

0.70 

0.71 

0.71 

Table  3 

.p- 

Values  of  4^2  ^or  various  values  of  n  and  f^ 
for  4*^*  =116  ft^  and  fj  =  0.046  sec"^. 


f 


u 

n 

0.46 

0.92 

4.6 

OO 

5.9 

.0034 

.0036 

.0038 

.0039 

5.7 

.0038 

.0042 

.0046 

.0048 

5.5 

.0044 

.0049 

.0057 

.0064 

5.3 

.0051 

.0060 

.0076 

.0102 

5.1 

.0060 

.0075 

.0105 

.0290 

4.9 

.0072 

.0097 

.0162 

00 

4.7 

.0087 

.0127 

.0261 

00 

4.5 

.0107 

.0167 

.0446 

00 
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Table  4. 

Values  of  1/Vtp^  for  various  values  of  n  and  f 

for  4'^*  *116  ft^  and  f,  »  0.046  sec’^. 
o  i 


f 

u 


u 

n 

0.46 

0.92 

4.6 

00 

5.9 

17.2 

16.7 

16.2 

16.1 

5.7 

16.2 

15.5 

14.8 

14.5 

5.5 

15.2 

14.2 

13.2 

12.5 

5.3 

14.1 

12.9 

11.5 

9.9 

5.1 

12.9 

11.5 

9.8 

5.9 

4.9 

11.8 

10.2 

7.9 

0 

4.7' 

10.7 

8.9 

6.2 

0 

4.5 

9.7 

7.7 

4.7 

0 

Table  5. 


Values  of  4'^^ 

/^o^2 

for  various  values 

n  and  f^ 

for  4<^* 
0 

=  116  ft^  and  f|  *  0.046  sec”^. 

f 

tt 

n 

0.46 

0.92 

4.6 

00 

5.9 

0.67 

0.63 

0.60 

0.59 

5.7 

0.63 

0.58 

0.53 

0.51 

5.5 

0.59 

0.52 

0.45 

0.40 

5.3 

0.54 

0.46 

0.37 

0.27 

5.1 

0.50 

0.40 

0.29 

0.10 

4.9 

0.46 

0.34 

0.21 

0 

4.7 

0.42 

0.29 

0.14 

0 

4.5 

0.38 

0.25 

0.10 

0 
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Di«cm«ion  of  tablet 

The  assumption  that  the  range  that  has  been  assumed  for  the 
values  of  f^  and  n  for  the  spectral  forms  for  S(k)  must  first  be  questioned. 
The  range  of  f^  is  probably  more  than  adequate.  However,  the  relation¬ 
ship  between  particle  motion  spectra  and  free  surface  spectra  is  not  known. 
Perhaps  n  can  be  even  smaller  than  4.5  and  still  yield  spectra  for  the 
free  surface  that  behave  over  a  given  range  like  where  p  is  greater 

than  5  as  wave  observations  seem  to  suggest. 

If  these  assumed  values  are  reasonable,  it  would  follow  that 
a  loop  in  a  realization,  z  =  z(x) ,  would  be  a  very  rare  event.  The  small¬ 
est  value  (except  the  zero)  in  Table  5  is  4.7  and  the  associated  proba- 
-5 

bility  is  less  than  10  If  loops  were  identified  with  breaking  waves,  the 
results  would  imply  that  waves  shorter  than  0.2  feet  and  n's  less  than 
5  produce  such  breakers  and  this  does  not  seem  at  all  reasonable. 

5 

Stated  another  way,  results  based  on  spectra  like  K/f 
depend  very  strongly  on  the  high  frequency  tail  of  the  spectrum  where  the 
waves  are  less  than  3  inches  long  and  on  the  fact  that  the  exponent  is 
exactly  5  and  not  5.1.  The  original  conditions  correspond  to  waves  with 
a  significant  height  of  42  feet  and  with  representative  lengths  of  1500 
to  2000  feet.  Breakers  in  such  a  sea  certainly  occur  and  surely  have 
dimensions  exceeding  3  inches.  Thus  although  such  loops  identify  un¬ 
realistic  spectra,  the  limitations  on  the  spectral  form  appear  to  be  due 
to  effects  that  occur  before  such  loops  even  have  a  chance  to  attempt  to 
form. 

The  tables  also  suggest  that  equation  (49)  will  be  approxi¬ 
mately  correct  as  to  general  shape,  and  that  the  higher  moments  are 
both  difficult  to  detect  theoretically  and  to  measure  in  practice. 

Equations  (45),  (46),  and  (52)  are  satisfactory  approximations,  but  the 
approximate  higher  moments  and  skewness  and  kurtosis  values  will  be  in 
error  by  large  percentages.  The  distribution  of  z  as  given  by  equation 
(44)  when  graphed  for  two  conditions,  one  for  which  f^  =  0.46 
and  n  =  5.9,  the  other  for  which  f^  =  4.6  and  n  =  4.5  can  hardly  be  dis  - 
tinguished  visually  from  the  normal  curve. 


Second  order  periodic  breaker  effect* 

The  Gerstner  wave  satiafies  the  nonlinear  Lagrangian  equations 
exactly.  However,  a  particle  never  moves  far  away  from  its  mean  position. 
A  strong  field  of  vorticity  exists.  The  limiting  steepness  of  a  Gerstner 
wave  is  given  by  ak  =  1.  For  irrotational  waves  Michell  [1893)  has  shown 

2&  2d. 

that  y-  =  0.142,  and  ak  =  1  corresponds  to  =  0.32  which  yields  a  wave 

much  steeper  than  is  possible  for  irrotational  motion.  It  was  mentioned 

above  that  the  Gerstner  wave  could  be  made  irrotational  to  second  order 

(with  third  order  complications)  by  adding  an  appropriate  drift  current. 

2 

With  6  equal  to  zero,  this  drift  current  becomes  a  uk. 

The  parametric  representation  for  the  free  surface  then  be¬ 


comes 

(62) 


X  =  a  -  a  sin^(ka  -  wt)  +  a^k  wt 
z  =  a  cos(ka  -  cot) 


for  a  periodic  wave,  which  can  be  transformed  by  solving  for  a  in  the 
second  equation  to  become 

(63)  cos"^(z/a)  -  ak  sin(cos*^(z/a) )  =  kx  -  cot-a^k^cot 


and  therefore  the  phase  speed  obtained  from  the  Lagrangian  equations 
solved  for  irrotational  motion  to  second  order  becomes 

(64)  C=t(l  +  aV). 

(l) 

The  phase  speed  obtained  by  solving  the  Eulerian  equations 
to  third  order  for  irrotational  motion  is  given  by 


The  Lagrangian  waves  to  second  order  are  thus  not  consistent 
with  Eulerian  waves  to  third  order.  The  second  order  periodic  irrotational 
Gerstner  waves  appear  to  be  traveling  as  if  the  linear  phase  speed  had 
been  simply  added  to  the  drift  current.  Third  order  periodic  irrotational 
Eulerian  waves  have  a  drift  current  that  is  the  same  as  the  Gerstner 
wave,  but  the  phase  speed  correction  is  only  half  of  the  phase  speed  cor¬ 
rection  for  Gerstner  waves.  The  discrepancy  can  be  explained  by  noting 
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that  there  probably  is  a  residual  vorticity  field  at  third  order  in  the  La- 
grangian  solution  that  may  account  for  the  difference. 

The  particle  velocity  at  the  crest  for  the  Gerstner  wave  is 

(66)  u(0.  0)  =  =  aa>  +  a^kw  . 

and  when  the  phase  speed  is  set  equal  to  the  particle  velocity,  the  limiting 
form 

(67)  ak  =  1 
is  still  obtained. 

With  ak  =  1.  the  second  order  phase  speed  and  the  particle 
velocities  at  the  crest  are  exactly  double  that  of  the  linear  theory,  but  it 
is  known  that  the  phase  speed  due  to  finite  height  cannot  exceed  about  1.2 
times  the  linear  phase  speed.  These  results  are  all  indicative  of  the  need 
for  different  wave  models  and  for  measurements  of  higher  order  effects 
in  waves. 

Construction  of  z  =  z(x)  for  the  second  order  equations 

From  recent  results  of  Tick  [1961],  it  may  be  possible  by  start¬ 
ing  with  equations  (28)  to  construct  the  vector  process  given  below  by  an 
appropriate  sequence  of  operations  on  the  linear  realizations  constructed 
above. 

x(o,  0,  0) 
z(a,  0,  0) 
x(o.  T,  0) 
z(a,  T,  0) 
x(a,  0,  6j) 
z(o,  0,  6j) 

x(o,  T, 6j) 

z(o,  T,  6j) 

x^(o,  0,  0) 


Xj(a,  0,6  j  ) 
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The  construction  may  even  turn  out  to  be  considerably  simpler 
than  that  employed  by  Tick  [1961]  because  only  low  frequencies  are  in¬ 
volved. 

In  the  above  process,  t  is  to  be  a  fraction  of  a  second  after 
t  =  0,  6j  is  to  be  a  foot  or  so  below  the  free  surface.  The  first  two 
pairs  will  yield  an  estimate  of  the  instantaneous  speed  of  a  crest.  If 
x^(a,  0,  0)  exceeds  this  value,  then  the  crest  could  be  considered  to  be  a 
breaking  crest.  The  other  terms  should  yield  some  information  on  the 
size  of  the  region  of  unstable  motion.  A  range  of  possible  spectral  forms 
can  be  used  to  construct  this  vector  process  and  some  idea  of  the  fraction 
of  breaking  waves  in  the  total  can  be  gained. 

However,  from  the  results  given  above  on  the  linear  random 
model  and  on  the  periodic  second  order  irrotational  model,  the  results 
obtained  from  such  a  construction  will  probably  not  be  quantitatively  cor¬ 
rect  as  they  may  yield  waves  that  are  too  high  physically  by  a  factor  of 
two,  even  for  irrotational  waves  to  second  order. 

Breakers  in  a  random  sea 

Whitecaps  and  breakers  in  deep  water  are  an  integral  part  of  the 
problem  of  wave  generation.  Many  photographs  of  waves  in  deep  water  show 
such  breaking  waves.  The  equations  governing  the  wave  motion  fail  locally 
and  a  breaker  is  produced.  The  water  in  the  breaking  part  of  a  wave  is 
certainly  governed  by  physical  laws  that  are  completely  different  from 
those  that  govern  the  wave  motion.  A  breaker  is  probably  produced  when 
the  water  particle  speed  exceeds  the  speed  of  the  crest,  but  other  effects 
such  as  the  wind  actually  blowing  off  the  crest  of  a  wave  before  this  re¬ 
quirement  is  met  may  also  play  an  important  part. 

Moreover,  the  turbulence  produced  by  a  breaking  wave,  governed 
by  eddy  viscosity  laws,  will  tend  to  damp  out  the  higher  frequency  linear 
components  in  the  spectrum  of  the  wind  generated  sea.  The  absence  of  the 
local  chop  in  the  wake  of  a  ship  illustrates  this  effect. 

Whenever  a  breaker  disorganizes  the  fluid  motions,  there  is 
a  mass  of  water  moving  forward  at  the  crest  of  the  waves.  This  effect 
may  contribute  to  both  the  drift  current  at  the  surface  and  to  the  growth 
of  the  lower  frequency  waves. 
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Alto  there  ii  a  wind  drift  current  in  a  local  sea  that  ia  nearly  in 
the  aame  direction  aa  the  traveling  wavea.  Thia  wind  drift  current  pro- 
ducea  a  field  of  vorticity  in  the  aame  aenae  aa  the  direction  of  propagation 
of  the  wavea  and  may  well  cause  breaking  to  occur  for  considerably  lower 
wave  heights  than  even  irrotational  theory  would  predict. 

The  work  of  Longuet -Higgins  [1953],  [I960]  has  shaken  the  con¬ 
cepts  of  irrotational  motion  for  gravity  waves  to  their  very  foundation  . 

His  results  predict  that,  even  in  the  absence  of  wind  drift  currents  and 
turbulent  effects  due  to  breakers,  the  mass  transport  velocity  in  deep 
water  may  be  considerably  stronger  than  irrotational  theory  predicts  if 
the  waves  have  been  running  long  enough. 

Moreover,  as  Longuet -Higgins  [1953]  has  pointed  out,  Duvriel- 
Jacotin  has  shown  that  any  mass  transport  velocity  can  be  used  as  the 
'  starting  point  and  a  wave  motion  can  be  superimposed  thereon. 

It  is  also  known  from  classical  periodic  wave  theory  that  a  sharp 
angle  at  the  crest  of  120*  is  another  way  to  define  the  limiting  height  of 
a  periodic  wave  of  the  form  z  =  z(x).  It  may  be  possible  if  the  other  con¬ 
cepts  discussed  above  fail,  due  to  the  fact  that  the  results  are  only  to 
second  order,  to  study  constructions  of  z  =  z(x)  as  a  random  process  to 
see  how  many  waves  approach  this  limiting  form. 

With  all  of  these  preliminary  remarks  it  would  seem  that 
equation  (29),  with  perhaps  still  an  arbitrary  second  order  vorticity  field 
added,  can  provide  some  insight  on  the  growth  of  a  wave  spectrum.  Con¬ 
sider  two  orbital  motion  spectra  with  the  same  variance,  one  with  con¬ 
tributions  with  relatively  higher  frequencies  than  the  other.  The  mass 
transport  at  the  surface  will  be  stronger  for  the  spectrum  with_the  richer 
high  frequency  content  but  at  the  same  time  the  speeds  of  the  crests  will 
be  less.  The  interaction  of  these  effects  could  conceivably  cause  the 
higher  frequency  waves  in  the  spectrum  as  they  are  overtaken  by  and 
ride  up  on  the  crests  of  the  longer  lower  frequency  waves  ( Longuet -Higgins 
and  Stewart  [I960])  to  break  and  dissipate  by  turbulent  action.  At  the 
same  time  the  low  frequency  components  can  continue  to  grow  and  even 
be  caused  to  grow  by  the  breaking  of  the  shorter  waves.  Considerations 

such  as  this  may  explain  various  wave  spectra  where  forms  like  K/w  , 

6  5  5  5 

K/u  ,  K/w  ,  and  K/u^  have  been  proposed  and  where  the  exponent 
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app«ars  to  depend  on  the  heights  of  the  waves  that  were  measured. 
Comparison  with  observations 

The  purpose  of  this  part  of  this  paper  is  to  present  the  results  of 
some  observations  that  suggest  that  the  various  models  proposed  above 
are  more  nearly  in  accord  with  actual  waves  than  models  based  on  the 
Eulerian  equations. 

Figure  4  is  a  copy  of  a  photograph  given  by  Roll  [1957]  for  waves 
a  few  centimeters  high  and  8  to  35  cm  long.  The  third  photograph  from  the 
top  shows  a  remarkable  similarity  to  figure  2.  For  various  wave  heights 
and  apparent  lengths,  either  a  rounded  crest  or  a  sharp  crest  occurs.  Note 
the  marked  departure  from  normality  that  slopes  and  curvatures  would 
have  if  they  were  evaluated  from  such  a  record. 

Figures  5,  6,  and  7  show  some  reproductions  of  seven  selected 
wave  profiles  obtained  by  stereo -photogrammetric  techniques  from  photo¬ 
graphs  taken  aboard  a  ship.  The  height  scale  is  five  times  that  of  the 
horizontal  scale.  The  profiles  show  the  presence  of  sharp  crests  and 
shallow  troughs  as  these  models  predict.  Note  also  that  the  profiles  are 
much  richer  in  high  wave  numbers  than  corresponding  time  histories  are 
rich  in  high  frequencies.  These  profiles  were  kindly  furnished  the  author 
by  Dr.  Norman  F.  Jasper  of  the  David  Taylor  Model  Basin.  Additional 
information  on  these  records  is  given  by  Brooks  and  Jasper  [1957]. 

Figure  8  shows  some  selected  time  histories  obtained  in  Buz¬ 
zards  Bay  by  Harlow  G.  Farmer  [Farmer  and  Ketchum  (in  press)].  The 
original  records  have  been  reversed,  so  that  time  increases  from  left 
to  right,  and  blackened  below  the  original  trace.  The  wave  sensing  sys  - 
tern  consists  of  a  very  fine  wire  (0.015  inch  diameter)  that  appears  to 
follow  the  water  elevation  more  accurately  than  previous  wave  poles  or 
wave  wires.  A  number  of  the  crests  are  very  sharply  pointed.  There 
is  evidence  that  a  breaker  occurred  in  the  top  trace  where  one  crest 
shows  a  nearly  vertical  rise  at  the  crest.  Since  these  records  are  much 
richer  in  high  frequency  content  (sharp  crests,  and  vertical  rises),  spectra 
computed  from  them  will  be  richer  at  high  frequencies  than  spectra  pre  - 
viously  reported  for  waves  in  this  general  range  of  heights  and  periods. 

The  figures  and  illustrations  in  Schulejkin  [I960]  also  suggest 
that  the  randomised  Lagrangian  model  corresponds  well  with  reality. 


Selected  wave 


The  three  modele  based  on  the  Lagrangian  equations  proposed 
in  this  paper  show  some  features  of  waves  that  are  observed  to  occur  and 
that  do  not  appear  to  be  present  in  models  based  on  the  EtUerian  equations. 
The  probability  structures  of  the  various  resulting  free  surfaces  are  not 
understood.  However,  in  the  simplest  model,  free  surfaces  can  be  con¬ 
structed  that  show  very  interesting  properties.  There  is  some  hope  that 
the  effects  of  breaking  waves  can  be  studied  by  means  of  these  models. 

It  is  also  evident  that  the  measurement  of  random  seas  needs  to  be  con¬ 
siderably  refined  so  as  to  detect  the  higher  order  effects  possible  in 
these  various  different  models. 

Acknowle  dgme  nts 

The  material  presented  in  this  technical  report  was  stimulated 
by  a  conversation  with  Dr.  Leo  J.  Tick  in  which  he  asked  whether  or  not 
random  wave  theory  could  be  approached  from  any  other  point  of  view 
than  that  presently  in  use.  If  breaking  waves  are  ever  to  be  studied  at 
all,  the  perturbation  procedure  in  the  Eulerian  system  of  equations 
does  not  seem  to  come  close  enough  soon  enough.  These  results  are 
a  start  toward  trying  to  understand  waves  from  another  point  of  view. 

The  help  of  Dr.  Tick  in  the  formulation  of  the  theoretical  developments 
is  also  appreciated. 

Stimulating  conversations  with  Dr.  G.  Neumann  with  reference  to 
eddy  viscosity  effects  and  visual  nonlinear  effects  of  breakers  were  most 
helpful. 


.38. 


i 


I 


Reference* 


Bieeel.  F.  [1952]:  Study  of  wave  propagation  in  water  of  gradually  va^rying 
de^h.  In  Gravity  Waves.  Natl.  Bur.  Standards  (U.  S. )  Circ.  521, 
pp.  221734. 

Brooks,  R.  L. ,  and  N.  F.  Jasper  [1957];  Statistics  on  wave  heights  and 
periods  for  the  North  Atlantic  Ocean.  DATMOBAS  Report  1091. 

Cartwright.  D.  E. ,  and  L.  J.  Rydill  [1956]:  The  rolling  and  pitching  of  a 
ship  at  sea.  Trans.  I.  N.  A.  (London),  vol.  98,  p.  100. 

Ehrenfeld,  Goodman,  et  al.  [1958]:  Theoretical  and  observed  results  for 
the  zero  and  ordinate  crossing  problems  of  stationary  Gaussian 
noise  with  application  to  pressure  records  of  ocean  waves. 

Tech.  Report  No.  1,  Bureau  of  Ships.  Contract  No.  Nobs  27018 
(1734F),  New  York  University,  College  of  Engineering,  Research 
Division. 

Farmer,  H.  G.  ,  and  D.  D.  Ketchum  [in  press]:  An  instrumentation  sys  • 
tern  for  wave  measurements,  recording  and  analysis.  Proc. 

7th  Conference  on  Coastal  En£ineering. 

Lamb,  H.  [1932]:  Hydrodynamic s ,  6th  Ed.  Cambridge  Univ.  Press. 

Lewis,  E.  V.  [1955]:  Ship  speeds  in  irregular  seas.  Trans.  SNAME, 
vol.  63. 

Longuet -Higgins,  M.  S.  [1953]:  Mass  transport  in  water  waves.  Phil. 

Trans.  Roy.  Soc.  of  London,  Series  A,  vol.  245,  pp.  535-581. 

Longuet -Higgins,  M.  S.  [1957]:  Statistical  properties  of  a  random  moving 
surface.  Phil.  Trans.  Roy.  Soc.  of  London,  Series  A,  vol.  250. 
pp.  157-171:  - 

Longuet -Higgins,  M.  S.  [I960]:  Mass  transport  in  the  boundary  layer 

at  a  free  oscillating  surface.  J.  Fluid  Mechanics,  8,  293-306. 

Longuet -Higgins,  M.  S. ,  and  R.  W.  Stewart  [I960]:  Changes  in  the 

form  of  short  gravity  waves  on  long  waves  and  tidal  currents. 

J.  Fluid  Mech.,  8,  565-583. 

Michell,  J.  H.  [1893]:  The  highest  wave  in  water.  Phil.  Mag.,  36, 

430-437. 

Phillips,  O.  M.  [1958]:  The  equilibrium  range  in  the  spectra  of 

wind -generated  waves.  J.  Fluid  Mechanics,  vol.  4,  part  4, 
pp.  426  -434. 


39- 


Phillips,  O.  M.  [1961]:  Th«  dynamics  of  random  finite  amplitude  gravity 

waves  (preprint).  NAS-NRC  Conference  on  Ocean  Wave  Spectra, 
May  1-4,  1961. 

Pierson,  W.  J.  [1952]:  A  unified  mathematical  theory  for  the  analysis, 
propagation  and  refraction  of  storm  generated  ocean  surface 
waves.  Parts  I  and  II.  Research  Division,  College  of  Engineer¬ 
ing,  New  York  University.  Prepared  for  the  Beach  Erosion 
Board,  Department  of  the  Army,  and  Office  of  Naval  Research, 
Department  of  the  Navy. 

Pierson,  W.  J.  [1955]:  Wind  Generated  Gravity  Waves.  Advances  in 
Geophysics,  vol.  2,  Academic  Press. 

Pierson,  W.  J. ,  and  P.  Fife  [1961]:  Some  nonlinear  properties  of  long 
crested  periodic  waves  with  lengths  near  2.44  centimeters. 

Journal  of  Geophysical  Research,  vol.  66,  no.  1. 

Roll,  H.  U.  [1957]:  OberflXchen-Wellen  des  Meeres.  Handbuch  der 
Physik,  (48),  Berlin,  Springer  Verlag. 

Schulejkin,  W.  W.  [I960]:  Theorie  der  Meereswellen.  Akad.  Verlas- 
Berlin. 


St.  Denis,  M.  ,  and  W,  J.  Pierson,  Jr.  [1953]:  On  the  motions  of  ships 

in  confused  seas.  Trans.  Society  of  Naval  Architects  and  Marine 
Engineers,  vol.  61,  pp.  Z80-357.  ^ - 

Stoker,  J.  J.  [1957]:  Water  Waves.  Interscience  Publishers,  New  York. 

Tick,  L.  J.  [1959]:  A  nonlinear  random  model  of  gravity  waves,  I. 

Jour.  Math,  and  Mech. ,  vol.  8,  no.  5. 

Tick,  L.  J.  [1961]:  Non-linear  probability  models  of  ocean  waves 
(preprint).  NAS  NRC  Conference  on  Ocean  Wave  Spectra, 

May  1-4,  1961. 


